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ABSTRACT 

We investigate the curvaton dynamics and the non-linearity parameters in 
curvaton model with potential slightly deviating from the quadratic form in 
detail. The non-linearity parameter g^L will show up due to the curvaton 
self-interaction. We also point out that the leading order of non-quadratic 
term in the curvaton potential can be negative, for example in the axion-type 
curvaton model. If a large positive g^L is detected, the axion-type curvaton 
model will be preferred. 



1 Introduction 



Inflation model [1] provides an elegant mechanism to solve many puzzles in the hot big 
bang model due to the quasi-exponential expansion of the universe. The quantum fluc- 
tuation generated during inflation naturally seeds the small anisotropies in the cosmic 
microwave background radiation and the formation of the large-scale structure. In gen- 
eral many scalar fields are present in the early universe. The dynamics of inflation is 
governed by inflaton(s), but the primordial curvature perturbation can be generated by 
inflaton(s) [2] or curvaton(s) [3], or both of them. However the feature of the primordial 
power spectrum caused by inflaton is similar to curvaton. 

The CMB non-Gaussianity [4] opens a new window to probe the physics of the early 
universe. A well-discussed ansatz of non-Gaussianity has the local shape. This kind of 
non-Gaussianity can be characterized by some non-linearity parameters: 

C(X) = C 9 (X) + jj/tf^W + ^LC 9 3 (X) + " " " , (1.1) 

where ( is the curvature perturbation in the uniform density slice, ( g denotes the Gaussian 
part of £. The non-linearity parameters J'nl and g^L parameterizes the non-Gaussianity 
from the irreducible 3-point and 4-point correlation functions respectively. 

In the case of single field inflation model j 1 ^ 1 ~ 0(n s — 1) [5], which is constrained 
by WMAP (n s = O.96OI0 013) [6] to be much less than unity. In some general inflation 
models [7] a large non-Gaussianity is possibly obtained, but the shape of non-Gaussianity 
is different from the local shape. As we know, curvaton model can easily generate a large 
local-type non-Gaussianity [3,8]. Different group reported different result of data analysis 
to f^f based on WMAP 3yr data: 

27 < fx? 1 < 147 (1.2) 

in [9] and 

- 70 < fjfi? < 91 (1.3) 
in [10] at 95% CL. WMAP group reported the 5yr result [6]: 

- 9 < fftf < 111 (95% CL). (1.4) 

Up to now a Gaussian distribution of the primordial curvature perturbation is still consis- 
tent with experiments. However, the lower bound on /^£ a/ in Eq. fll.4l) is greatly reduced 
from WMAP3 to WMAP5. If a large positive /^£ a/ is confirmed by more years inte- 
gration of WMAP, or the Planck mission, it can help us to distinguish curvaton model 
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from inflation model. Recently many issues related to curvaton were widely discussed 



Curvaton is a light scalar field with weak self-interaction. Even though its energy den- 
sity is subdominant during inflation, it can make the main contribution to the primordial 
curvature perturbation. In curvaton model, J'nl is inversely proportional to the fraction 
of curvaton energy density in the energy budget at the epoch of curvaton decay, and the 
value of g^L is mainly contributed by the curvaton self- interaction. In the literatures, the 
interaction term of curvaton is assumed to be neglected and then g^L is small. Now the 
cosmological observations become more and more accurate, and even g^L can be detected 
in the near future if it is not too small. So it is worth working out how g^L is related to 
the curvaton self-interaction. 

In [13,14] the authors only considered the case with small positive interaction term and 
they concluded that g^L ~ C(— 10 4 ) — 0{— 10 5 ) for the reasonable cases. In this paper 
we will see that it is also reasonable to consider the case with small negative interaction 
term for curvaton model and discuss the dynamics of curvaton more carefully. 

This paper is organized as follows. In Sec. 2, we discuss some possible forms of the 
curvaton potential. In Sec. 3, the dynamics of curvaton with potential slightly deviating 
from quadratic form is studied in detail. In Sec. 4, we treat the non-quadratic potential 
as a perturbation, and analytically calculate f^L an d gNL to the leading order. How to 
determine the non-quadratic term from experiments will be investigated in Sec. 5. At the 
end, we give a summary in Sec. 6. 

2 Estimator of small deviation from the quadratic 
potential 

We consider the following form of the curvaton potential, 



The mass term is assumed to be dominant, and the size of the interaction term compared 
with the mass term is measured by 



where a* is the value of curvaton at Hubble exit. In curvaton scenario the curvaton mass 
m is assumed to be much smaller than the Hubble parameter if* during inflation, and 
then the value of curvaton at the end of inflation is also o~*. 
1 Here the value of s is half of that in Ref. [13, 14]. 



in [11-13]. 




(2.1) 




(2.2) 
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Usually the stability of the system requires the coupling constant A to be positive. 
Otherwise the potential does not have a lower bound when the vacuum expectation value 
of the field goes into infinity. The authors in [13, 14] only focused on the case of s > 0. 
However we cannot expect the potential is still reliable when a is very large. So we can 
expect that s can take a negative value. 

On the other hand, to keep the mass of curvaton small enough, some symmetries 
are called for to prevent the large quantum correction to its mass. It seems natural to 
invoke supersymmetry. However supersymmetry must be broken down during inflation 
and the mass square of each scalar field generically receives a correction of order H%. An 
alternative candidate of curvaton is the pseudo-Nambu-Goldstone boson [15] - axion. 
The potential of axion is given by 



where / is called the axion decay constant. The smallness of the axion mass is protected 
by the shift symmetry a — > a + 5. In string compactifications, axion fields are popular, 
and even when all other moduli are stabilized, the axion potential remains rather flat as a 
consequence of nonrenormalization theorems. In the inflation driven by axion, the decay 
constant is required to be larger than the Planck scale. However such a large value of / 
cannot be achieved in string theory [16,17]. But the axion is a nice candidate of curvaton. 
For a small axion displacement a < f, the potential is expanded as 



The mass and the coupling of axion are respectively m = M 2 /f and A = -M 4 /(24/ 4 ). 
In this case the parameter s is 

which is negative. Since a* < /, —1/24 < s < 0. 

Because the curvaton mass is small compared to the Hubble parameter if* during 
inflation, the Compton wavelength of curvaton is larger than the background curvature 
radius H' 1 , and the gravitational effects may play a crucial role on the behavior of 
curvaton. In [18] the authors explicitly showed that the quantum fluctuation of a light 
scalar field a with mass m in de Sitter space gives it a non-zero expectation value of a 2 
as follows 




(2.3) 




(2.4) 




(2.6) 
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So the typical value of curvaton at the Hubble exit is given by 



(2.7) 



In 2 m 



and then 



For the axion-type curvaton, we have 

1 Hi 1 Hi 




(21 



(2.9) 



(87r) 2 m 2 / 2 (8vr) 2 M 4 ' 
Here we assume er* < f. Let's consider a concrete example of axion in [16] where a 
five-dimensional U(l) gauge theory is compactified on a circle of radius R. A non-local 
potential as a function of the gauge invariant Wilson loop 

e ie = e ifA 5 dx5 ( 21 y 



£= o 2^, (dO) 2 -V(9), (2.11) 



will be generated in presence of charged fields in the bulk. The effective Lagrangian for 9 
takes the form 

1 

2gl(2~nRy' 

where the four-dimensional gauge coupling g± is related to that in five dimensions by 
g\ = fi , 5/(27ri?), and the potential is roughly given by 

V{e) ~lt^ {1 - C0S6) ' (2 ' 12) 
and c ~ (9(1) depends on the number of charged fields in the bulk. In this case, / = 2 J gi R i 
M 4 = Y^-^and the parameter s becomes 

s = -^(H*R) A . (2.13) 

In order for the validity of the effective description of the four- dimensional scenario, the 
size of extra dimension is required to be much smaller than H~ x . If s can be determined 
by experiments, we can roughly know the ratio between the size of extra dimension and 
the Hubble radius during inflation in this model. 



3 The dynamics of curvaton after inflation 

After inflation radiation dominates our Universe and the Hubble parameter drops with 
time as H = l/(2t). The curvaton equation of motion after inflation is 

a + —a = -m 2 a ( 1 + n\(-) n - 2 ) . (3.1) 
2t \ m / 
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The correction in the equation of motion is small if \s\ <C 1/n. We introduce a dimen- 
sionless time coordinate x = mt and the curvaton equation of motion becomes 

// 3 , 



2x 



-a 



-a 



(l + n\(-: 
V m 



n-2 



(3.2) 



here the prime denotes the derivative with respect to x. The time when inflation ends is 
roughly given by t en d = 1/(2H*) and thus the initial time coordinate Xi n i = rnt en d — 0. 
Assume x = x Q = mt Q at the time of curvaton starting to oscillate. Perturbatively solving 
Eq. (13.21) with initial conditions a^i = a* and a' ini = at x = 0, we find the value of 
curvaton at x n is 



(To = cr* 



w(x ) + 77,^(77,, x )X( 



CT 



* \n-2 



m 



0{s 7 



where 



w[x ) 
gin, x ) 



2ir(5/4)x 3 Ji(x ), 

n — 5 

tt2— r(5/4) ? 



\n-l 4 



X 



Ji(x ) I J" 1 (x 



Ji(x)x t dx 



(3.3) 



(3.4) 



(3.5) 



'0 4 Jo 
and Jx/4 and Y1/4 are respectively the 1/4 order Bessel J function and Bessel Y function. 
The curvaton starts to oscillate when the Hubble parameter drops below the same order 
of magnitude of the curvaton mass m. Usually we have H = m and then x — ~. 
In [14] the authors assumed x Q — 1. The values of w(x ) and g(n,x Q ) for different x Q are 



showed in Fig. 1. For x Q 



2 ■ 



w(x 



0.95 and g{2,x Q ) ~ -0.049, g(±,x ) 



-0.047, 



w(x o 




9(n,x, 




Figure 1: The values of w(x Q ) and g(n,x ) for different x a . 

g(6,x ) ~ -0.046. For x c = 1, w(x Q ) ~ 0.81 and ^(2,x c ) ~ -0.179, ^(4,x G ) ~ -0.16, 
^f(6, x D ) ~ —0.145. We see that w(x ) and g(n, x Q ) depend on the choice of the time when 
the curvaton starts to oscillate. 
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4 Non-linearity parameters in curvaton model 



The primordial curvature perturbation was calculated in [3,8]. The primordial power 
spectrum, bispectrum and trispectrum are defined by 

(C(k 1 )C(k 2 )) = (2n) 3 V c (k 1 )5 3 (k 1 + k 2 ), (4.6) 
(C(k 1 )C(k 2 )C(k 3 )) = (2nfB c (k 1 , k 2 , fc3)<f 3 (k! + k 2 + k 3 ), (4.7) 
(C(k 1 )C(k 2 )C(k 3 )C(k 4 )) = (27i) 3 T c (k 1 , k 2 , k 3 , fc 4 )* 3 (ki + k 2 + k 3 + k 4 ). (4.8) 

The bispectrum and trispectrum are respectively related to the power spectrum by 
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B c {h, k 2 , h) = -fNLlV^k^ih) + 2 perms], 

T c (ki,k 2 ,k 3 ,k 4 ) = TNLlV^ky^V^k-^V^k^) + 11 perms] 
54 

+ ^gNL[V^k 2 )V c (h)V^h) + 3 perms]. 



(4.9) 
(4.10) 



(4.11) 



Here the non-linearity parameter tnl depends on f^L- 

r - 36 f 2 

~ 25 

In curvaton model, the amplitude of the primordial power spectrum and the non-linearity 
parameters can be written as 



9vr 
5 



1 f 2 2 



4/ 
25 

9NL ~U 



h) 



D 



5 
3 



5/d 
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+ 3h) -±(l + h) + hl- 9h) + 10f D + 3f D 



where 



h 



3U 



<7,D 



4- a 



< 2 ' < 3 



According to Eq. (13.31) . the parameters q, h and h take the form 
w{x ) + n{n — l)g(n, x Q )s 



Q 
h 
h 



w(x Q ) + ng(n, x Q )s 
w(x ) + ng(n,x )s 
{w(x ) + n{n — l)g(n, x )s) 2 

(w(x ) + ng(n,x )s) 2 
{w(x ) + n{n — l)g(n, x )s) 3 



n(n - l)(n - 2)g(n,x )s, 

n(n — l)(n — 2){n — 3)g(n, x )s. 



(4.12) 
(4.13) 
(4.14) 

(4.15) 
(4.16) 

(4.17) 
(4.18) 
(4.19) 



We see that both h and h are proportional to the size of the non-quadractic term s. If 
the curvaton potential deviates from the exactly quadractic form, the value of Jnl can 
be small even when fjj <C 1 if 1 is cancelled by h, but g^L is still very large. 

Here we also want to stress that the choice of the time when curvaton starts to oscillate 
has much effect on the quantitative results. For n = 4, we illustrate the quantitative 
results of q, h and h for different choice of x Q in Fig. 2. 




Figure 2: The quantitative results of q, h and h for different x Q . Here we take n = 4. 

In [13] the authors only consider the case of s > and they concluded that /tvl can be 
small, but g^L will be 0{— 10 4 ) — 0{— 10 5 ) for the reasonable cases. In Sec. 2 we argued 
that it is reasonable to extend the discussions into the case of s < 0. For s < 0, both 
h and h are positive, /nl always gets a large positive value if fo 1, and g^L ~ Jnl 
which is roughly the same order of magnitude as t nl . See the numerical results of f^L 
and g NL in Fig. 3. 
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Figure 3: The non-linearity parameters f^L and g^L for different n, s and x Q are plotted. 
Here we take fu = 0.01. 

5 A proposal to measure the non-quadratic term 

In [13] the authors pointed out that any deviation from the relation 

9nl ^ -y/Wi (5.1) 

will indicate that the curvaton potential does not take the exactly quadratic form. Here 
we will take a closer look at the measurement of the curvaton self-interaction term. 

From Eq.(jH3]) and (Ol) . we find 

10 , 25H3/i 25, , , , , 

9NL + ~f NL = YC-f^~ + T08 ( " 23 ~ 9k + 8f ° + 6fD) - (5 ' 2) 

Here we only focus on the case with large absolute values of /nl and g^L, because the 
accuracy of experiments is limited. Or equivalently, we have /d <C 1. If \s\ <C 1, we have 
Jnl and 



10 2 

9nl + -z-fNL - o / \ fNL^ 2 {n - l)(n - 2)g(n,x )s. (5.3) 

3 6W{X ) 
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For convenience, we define a new quantity a(s) as follows 

{ } ~ 2p NL n*(n-l)(n-2)g(n,x o y 1 " ] 

If |s| <C 1, we can expect a(s) = s. Since there are three free parameters (fo, n, s) from 
the theoretical viewpoint, but only two observed quantities (Jnl, 9nl), we cannot model- 
independently determine the value of s from experiments. For an instance, we consider 
the interaction term takes the form Act 4 and show a(s) as a function of s in Fig. 4. Even 




or(s) 



Figure 4: This figure illustrates a(s) as a function of s. In fact, for the large value of Jnl, 
it is insensitive to f^L- 

though we set f^L = 50, one can check that a(s) is quite insensitive to f^L- In this case 
a(s) ~ s if \s\ < 0.01 and this result is roughly independent on the choice of x Q . 

Before the end of this section, we also want to estimate the order of magnitude of the 
dimensionless coupling coupling A in the case of n = 4. The typical value of curvaton 
during inflation is ct* ~ H^/m, and then s ~ \(H*/m) 4 and the amplitude of the primor- 
dial power spectrum is ~ /|,m 2 /if 2 . For large /atl and small s, we have /nl ~ V/d- 
Therefore m/H* ~ ^/P^/nl and A ~ s(y/P^fNL) i - Taking into account WMAP nor- 
malization P c = 2.457 x 10~ 9 [6], A ~ 10~ 18 f^ L s. For f NL ~ 10 2 and s ~ ±10~ 2 , 
Qnl ~ _ IOs/atl ~ =Fl0 3 and A ~ ±10~ 12 . Such a weak interaction can be possibly 
detected by the cosmological experiments in the near future! 
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6 Conclusions 



The fluctuations of curvaton field evolve non-linearly on the superhorizon scales if the 
curvaton potential deviates from the exactly quadratic form. In this paper, we suggest 
that the leading order of the non-quadratic term for curvaton field can be negative, for 
example in the axion-type curvaton model. We also investigate the curvaton dynamics 
and the non-linearity parameters in curvaton model with non-quadratic curvaton potential 
in detail. For a positive non-quadratic term, /jvl can be very small even when /o<l. 
But it does not happen in the case with a negative non-quadratic term. We also see that 
the second order non-linearity parameter g^i measures the size of the self- interaction of 
curvaton field. If g^L is positive, the leading order of the self-interaction term should 
be negative, and the axion-type curvaton model is preferred. The next generation of 
experiments such as Planck will improve the accuracy to A/tvl ~ 6 and At/vl ~ 560 [19]. 
We hope the non-linearity parameters will be detected soon. 

In this paper, the self-interaction term is treated as a perturbation. The case in 
which the self-interaction term dominates the curvaton potential during inflation will 
be discussed in [20]. On the other hand, multiplicity of curvaton fields is expected in 
the fundamental theories going beyond the standard model. In particular, axions are 
typically present in large numbers in string compactifications. In [12], these axion fields are 
suggested to be taken as curvatons. It is worth investigating the non-linearity parameters 
in N-vaton model [12] carefully in the future. 
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